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Abstract 

The renormahzation issue of the reparameterization invariance in heavy quark effective theory 
and NRQCD is investigated. I argue that the renormahzation of the transformation of the heavy 
quark field under the variation of the velocity parameter V is attributed to the renormalization of 
the small component field in the proposed transformation. I show that the matching condition for 
determining the renormalized small component field can be obtained by imposing an infinitesimal 
transformation of V on the relations between the Green's functions in QCD full theory and those 
in the effective theory. As an application, I determine the renormalized transformation to order 
l/m? using the matching condition. The obtained result is in disagreement with that determined 
by indirect method. 
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I. INTRODUCTION 

Heavy quark effective theory (HQET)^ and nonrelativistic QCD (NRQCD)P| are pow- 
erful tools in dealing with dynamics of heavy-light and heavy-heavy systems, respectively. 
In those systems, the off-shell momentum of the heavy quark is much smaller than its mass. 
The effective theories are designed to reproduce the results of the QCD full theory at the 
low energy scale in a simpler way by integrating out the effects at the energy scale of the 
heavy quark mass. In the past decade both effective theories and their applications have 
been intensively studied. 

One interesting theoretical issue in those effective theories is the reparameterization in- 
variance (RPI). It arises from the fact that the effective theory explicitly depends on the 
four velocity parameter V. In constructing the effective lagrangian, one needs to divide the 
heavy quark momentum P into a large part and a small one as P = mV + k, where m is the 
heavy quark mass and is a small residual momentum, one also needs to decompose the 
Dirac 4-fermion field as large and small two component fields in respect of V and uses the 
large one to describe the heavy quark or antiquark. These procedures lead to the effective 
lagrangian being V^-dependent. The choice of V which satisfies V^^ = 1 is not unique. But 
the physical prediction should be unchanged against the variation of the velocity parameter 
V. This is the RPI. It is required by the consistence of the effective theory and also conducts 
interesting applications. It was first proposed in HQET. However, the same invariance also 
holds in NRQCD effective theory. 

To implement RPI in the effective theory, it is essential to find out an appropriate trans- 
formation of the heavy quark field under the variation of V, which is found to be quite 

n 

nontrivial. It was first studied by Manohar and Luke^] in HQET. They used the Lorentz 
boost of spinor field as the transformation of the heavy quark field from finite velocity 
V —>■ v. Their transformation suffers from operator ordering ambiguities when it is ex- 
panded to a higher order of 1/m while Manohar discussed its expansion at higher order. 
Chen I proposed an infinitesimal transformation of the heavy quark field under the velocity 
variation from V ^ V + AV. Chen's transformation keeps tree level effective lagrangian in- 
variant to all orders of 1/m. Finkemeier, Georgi, and Mclrvinj^ showed that to order 
the effective lagrangian constrainted by Manohar and Luke's transformation and Chen's 
transformation may be related by a field redefinition. 
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Chen's transformation can be expanded as inverse power series of the heavy quark mass. 
Each term contains the product of some covariant derivatives and the heavy quark field. 
Since this kind of expansion changes the ultraviolet behavior of the original transformation, 
beyond the next leading order, it turns out that it needs to be renormalized. The renor- 
malization of it is different from that of the effective lagrangian in which each term is the 
bilinear function of the heavy quark fields. For that case, we have appropriate matching 
conditions to determine the renormalized transformation. 

Kilian and Ohl[?] proposed a renormalized transformation. The form is exactly the same 
with Chen's transformation except that the covariant derivative in Chen's transformation 
is substituted by another operator which they called as the general covariant derivative. 
Sundrum^l discussed this issue using the auxiliary field method and obtained a result similar 
with Kilian and Ohl's. The results presented in these papers are formal. They did not showed 
how to determine the general covariant derivative by some specific matching conditions. 
Actually, in literatures no specific calculation for determining the transformation has been 
done with this method. The only calculation to determine the renormalized transformation 
is given by Balzereit in an unpublished paper P]. He first calculated the effective lagrangian 
to order at one loop level in the leading logarithmic approximation. By requiring 

the effective lagrangian invariant, he then determined the renormalized transformation to 
order indirectly. There are some drawbacks in this kind of calculations. First, they 

are quite complicated since the determination of higher order effective lagrangian usually is 
hard work. Second, it makes the RPI less practical application. An interesting application 
of the RPI is that once we know the transformation, we can use it to constraint the higher 
order effective lagrangian and makes the calculations simpler . Balzereit did it in an 

inverse order so RPI may not be used to constraint the effective lagrangian in his method. 

In this paper, I study this issue in an alternative way, with emphasizing on determining 
the renormalized transformation using the matching conditions. It is essential in study- 
ing the renormalization issue of RPI. With it, one is able to determine the renormalized 
transformation to any desired order both in as and in 1/m expansion. 

To derive the matching condition, it is instructive to recall the renormalization procedure 
of the effective lagrangian. There are simpler relations between the Green's function in the 
unexpanded effective theory with nonlocal form. These relations ensure that the effective 
theory and the full QCD reproduce the same physical predictions. The expanded effective 
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theory changes the ultraviolet behavior of the unexpanded effective theory. In order to 
reproduce the same result of the full QCD, those relations are required to be satisfied and one 
has to add some local operators in the effective lagrangian, i.e., renormalizing the effective 
lagrangian. Therefore, those relations between the Green's functions can be used as the 
matching conditions 1] to determine the effective lagrangian. On the other hand, there is no 
direct relations between the IPI vertexes in the effective theory and those in the full QCD. 
Thus one may not obtain the matching conditions from the IPI vertex. 

As we will see later on, the renormalization of the transformation of the heavy quark field 
under the variation of the velocity parameter V is attributed to the renormalization of the 
small component field in the proposed transformation. Therefore, if we can find some rela- 
tions between certain correlation functions in the full QCD and those in the effective theory, 
which involve the small component field, it then can be used as the matching conditions 
to determine the renormalized transformation. It is found that imposing an infinitesimal 
transformation on both sides of the relations between the Green's functions in QCD full 
theory and those in the effective theory, new relations in which the small component field is 
involved on the effective theory side are obtained. These new relations are nothing but the 
matching conditions to determining the renormalized small component field. They provide 
a systematic way to determine the renormalized transformation. As a specific example, I 
use these matching conditions to determine the renormalized transformation to order 
at one-loop level. The obtained result is in disagreement with that obtained by Balzereitj^. 
Since the same lagrangian (MRR lagrangian) are used, the disagreement between these two 
different results cannot be accounted by a field redefinition. I then show that the renormal- 
ized transformation determined by these matching conditions can be written in the form of 
Chen's transformation with the covariant derivative substituted by an operator which may 
be called as a general covariant derivative. Thus renormalized transformation determined by 
.atcM„. condMo. p.e.,.ed . tM= pape. i. con..en. with KUian a.d OU, ..ult^ 

The remainder of the paper is organized as follows. In section II, after a brief review 
on the tree level transformation, I argue that the renormalized transformation of the heavy 
quark field is attributed to the renormalization of the small component field. I then show 
that the matching condition for determining it can be obtained by imposing an infinitesimal 
transformation on the general relations between the Green's functions in QCD full theory 
and those in the effective theory. As an example, in section III, I determine the renormalized 
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transformation to order by matching the two-point and the three-point functions. In 

section IV, I show that previous resuhs can be more clearly understood in an alternative 
way to construct the effective lagrangian, where a four-component effective lagrangian con- 
structed first, followed by its reduced into the effective lagrangian in a two-component field. 
I then show that the remormalized small component field determined by the matching con- 
ditions are consistent with Kilian and Ohl's result.}^. Then I shown that the remormalized 
effective lagrangian is reparameterization invariant under the renormalized transformation. 
Section V contributes to the conclusion. Finally, in appendix A, I derive the general relation 
between the Green's functions using the generating functional method. 



II. RENORMALIZED TRANSFORMATION OF THE HEAVY QUARK FIELD 

In the heavy quark effective theory, the heavy quark is described by a two component 
field while in QCD full theory it is described by the Dirac four-component field. Thus, 
to construct the effective theory, one needs first to decompose the Dirac four- component 
field into two-component field. A simpler way to realize this decomposition introduced by 
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is: 

hY±{x) = exp{imV-x) P±\l'(x), (1) 



where 

^ _ 1± f 



(2) 



2 ' 

are the projection operators. The phase factor just removes the large part of the heavy 
quark momentum when it is written as p = mV + k. This definition of the field is used 
by most people [l| in literatures. Nevertheless, it is not unique. Different definitions lead to 
different forms of the effective lagrangian. However, they can be related to each other by 
field redefinition and produce the same physical predictionsj;^] . Overall in this paper, we use 
the definition of 



A. the tree level transformation 



With the definition (^J, the effective lagrangian reads 



= hv+{x)zD-Vhv+ix) -hv+ix) P J Pv-Phv+{x), (3) 

2m + iJJ ■ V 
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where D'^ = —iQsA^ is the covariant derivative. This is the nonlocal form of the effective 
theory. 

Obviously, this effective lagrangian depends on the velocity parameter V . The choice of 
this parameter is not unique. The RPI implies that physical predictions by the effective 
theory are independent of the choice of V. In Ref. |5i|, it was shown that the effective 
lagrangian Q is invariant under an infinitesimal transformation V ^ V + AV 

Ahv+{x) = [hv+{x) + hv-{x)) . (4) 

with the hy_{x) being the small component field and given by: 

AV is constrained by AV ■ V = due to V^^ = 1. 

Both the effective lagrangian (jS)) and the transformation given by (0)) and (0) can be 
expanded as power series of 1/m. The RPI is then valid order by order in 1/m. It implies 
that the tree level transformation makes the tree level effective lagrangian valid at any 
order of 1/m. The cancellation of the lagrangian shift AL at each order is quite nontrivial. 
Without expansion, the effective theory with the lagrangian is equivalent to that of the 
full theory in the sense that they produce the same S*— matrix elements. 



B. matching condition for renormalizing the transformation 

With expansion in terms of 1/m, the ultraviolet behavior of the theory is changed. To 
compensate this change, both the effective lagrangian and the transformation receive renor- 
malization. The renormalization procedure of the effective lagrangian is well-known while 
it is not evident how to renormalize the transformation. The difficulty arises from the fact 
that each term contains the product of some covariant derivatives and the heavy quark 
field which is complicated and it is just the linear function of the heavy quark field but not 
bi-linear functions as those in the effective lagrangian. 

Nevertheless, to gain the answer, it is instructive to recall the renormalization procedure 
of the effective lagrangian. For the heavy quark field defined in (^, there are simpler relations 
between Green's functions in the effective theory with the nonlocal form (jH)) and those in 
QCD full theory. For the expanded effective lagrangian, one has to add some local operators 
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to the expanded effective lagrangian so that those relations can be satisfied. Consequently, 
the renormalized effective lagrangian can reproduce the results of full QCD. Therefore, those 
relations between the Green's functions are just the matching conditions to determine the 
effective lagrangian. 

In order to renormalize the expanded transformation, the key point is to gain appropriate 
matching conditions to determine the transformation. As we will see, the renormalization 
of the transformation of the heavy quark field under the variation of the velocity parameter 
V is attributed to the renormalization of the small component field in the proposed trans- 
formation (j3]). Therefore, if we can establish some relations between the Green's functions 
in QCD full theory and certain correlation functions which involve the hy_{x) field in the 
effective theory, they can be used as the matching conditions to renormalize the hv-{x) 
and hence the transformation. It is found that these relations can be obtained by impos- 
ing an infinitesimal transformation of V on both sides of the relations between the Green's 
functions in QCD full theory and those in the effective theory. 

Now let's first look at the relations between the Green's functions in the full QCD and 
those in the effective theory. 

Denote the Green's function in the full theory by G{x,y;B) and that in the effective 
theory by Gv{x,y; B), respectively, where B is an arbitrary background field. They are 
defined by 

G{x,y;B) = ( | T ^(y) | )^ , (6) 

and 

Gv{x,y;B) = {0\T hv+{x)hv+{y)\Of , (7) 

Any interaction vertex with gluon can then be obtained by functional differentiating over 
B{x). When the quark field hv+{x) is related to the field in the full theory by Eq.((T)), the 
Green's functions satisfy the following simple relation (A derivation of this relation using 
generating functional method is given in Appendix A): 

Gv{x,y;B) = P+ 5) P+ , (8) 

here = means that we omit the phase factor exp [imV ■ {x — y)] and the renormalization 
constant Z{m, as{m)) which arises from the renormalization of the heavy quark field in the 
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full theory and that in the effective theory. Both sides are valid to all orders of 1/m and 
Us- This relation ensures that the S*— matrix elements in the effective theory are identical 
to those in the full theory. The local effective theory is gained by expanding the nonlocal 
effective lagrangian density Q. This expansion changes the ultraviolet behavior of the 
nonlocal effective theory. To reproduce the result of the nonlocal effective theory, one needs 
to add local operators to the expanded effective lagrangian. The relation (jSJ is required to 
be satisfied as matching conditions to determine the coefficients of those local operators. 
Thus the relations between the Green's functions in the full theory and those in the effective 
theory can be used as the matching conditions to determine the effective lagrangian. 

Below starting from these relations we derive the matching conditions to determine the 
renormalized transformation of the heavy quark field. 

Obviously, these relations are valid for arbitrary V . It allows us impose an infinitesimal 
transformation V + AV on both sides. It follows that: 

AGv{x,y-B) = -^G{x,y-B)V+ + Y+G{x,y-B)^ . (9) 

Again the symbol = means we omit the phase factor and a term arises from its infinitesimal 
shift which is trivial under the transformation. 

Given the definitions of the Green's functions in Eqs.©,©, we have the following unique 
solution of 0: 

Ait 

(o|r(A/^v+(x)/iv+(i/))|o)^ = ^(o|rv[/(x)t(i/)|o)^p+. (lo) 

Eq. (fTUj) implies that A/iy_|_(x) is proportional to .A^. Thus we can generally write it as 
the following form 

A/Iy+(X) = ^ (P+ h'y_,{x) + P_ h'y_{x) ) . (11) 

Substituting it into (fTUj). Eq. (fTUj) is then decomposed into two equations by projection 
operator: 

{Q\T{h'y^{x)hy^{y))\Q)^ = P+(0|TvI/(x)^(y)|0)^P+ , (12) 

{o\Tih'y_ix)hvAy))\o)'' = P-(o|rvi/(x)t(i/)|o)^p+. (13) 

Eq. (fT^ is nothing but Eq. (jH)) if hy_^^{x) is identical to hv+{x). Eq. (fT^ is a new one which 
can be regarded as the definition of the h'y_{x). Some comments can be given about this 



equation. The right hand side is still the Green's functions in the full theory while the left 
hand side is the Green's functions in the effective theory with insertion of local operators 
at point X. Again both sides are valid at any desired order of 1/m and with all possible 
interaction with gluons. 

For the effective lagrangian at tree level, this equation is satisfied if A/iy+(x) is given by 
transformation (jH) with hv-{x) given by the expanded expression At loop level, it is 
also valid at any specific loop momentum which is smaller than the quark mass. However, 
in calculating the whole loop momentum integration, some differences arise and the (|T^ can 
be regarded as the matching condition to determine the renormalized hv-{x) field. 

Eq. (jlip together with (|12p and (|13|) imply that the renormalized transformation keeps 
the same form as the tree-level transformation (^. But only the small component field needs 
to be renormalized. 

Now let's illustrate how p3|) determine the renormalized hy_{x) field as the matching 
conditions. For simplicity, we limit our arguments in the hard cutoff regularization. Similar 
arguments are applicable to dimensional regularization. 

Suppose one takes different hard cut-off regularization energy scales Ag and A/ in effective 
theory and in full theory, respectively. The A f should be much larger than the heavy quark 
mass m for including both quark and antiquark contributions. The Ae should be much 
smaller than m for the validation of the 1/m expansion. So they satisfy a hierarchy relation 
Ae ^ m ^ A/. Thus in calculating loop diagrams the integration bounds are different on 
both sides ( the left hand side is integrated out from to Ag while the right hand side is 
integrated from to Aj ). It leads to that the tree level expression of hy-{x) given by the 
expansion of (0) makes relation (|T^ no longer valid. To compensate the differences, one has 
to add the contributions of the loop momentum integrals from Ag to Aj to the left hand side. 
Those contributions can be expressed as the insertions of the local operators on the heavy 
quark line. They may be at point x or not. Those local operators at point x can be absorbed 
into the redefinition of the h'y_{x) while those local operators not shrunk at the point x 
correspond to the insertions of the higher dimensional operators in the effective lagrangian. 
Therefore, Eq. (fTT|) is just the matching conditions for determining the renormalized small 
component field hy^{x). Eq. ()13p allows one determine the renormalized hy_{x) to any 
order of 1/m and a^. 
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III. RENORMALIZED TRANSFORMATION TO ORDER l/m^ 



As a specific example, in this section, we determine the renormahzed effective lagrangian 
and the hv-{x) field up to the next leading order corrections of 1/m using the matching con- 
ditions (jHl) and (jl3p . respectively. Here we use the dimensional regularization and Feynman 
gauge in all specific calculations. 

Up to next leading order corrections of 1/m, the most general form of the renormalized 
effective lagrangian can be expressed as: 

Ly{x) = Z hv+{x)iD -V hv+{x) + -^hv+ix)'D'^ hY+{x) 



ZZ - zz 
+ ^ hv+{^) {^D ■ Vf hy+{x) + hy+{x) a^'^G.^ix) hy+{x) , (14) 

^^i.u ^ g^j^j most general form of the renormalized hv-{x) field can be written 



hy.{x) = F^^tP + ^D.V P + ^ PD.V)hy^{x), (15) 



as: 

2m Am? 4m^ 

where Z, Zg, Zm, do{fi), di{fi), (i2(yu) are the short distance coefficients to be determined. 

In the following matching procedures, all the short-distance coefficients are assumed to be 
calculated to all orders in a^, which makes the matching procedure are applicable for higher 
order calculations. But in this paper we only present one-loop result in the final expression. 

Since the off-shell momenta of the heavy quark and the momenta of the gluons can be 
treated as much smaller than the quark mass, in the matching procedure, the integrand can 
be expanded as power series of these momenta over the heavy quark mass. It leads to that 
the remainder part of the loop momentum integrals no longer depends on the heavy quark 
mass. We are free to choose the infrared regulator since the infrared divergences cancel on 
both sides. If we use a limitation order in which the external off-shell momenta of heavy 
quark and the momenta of the gluons go to zero first, followed hj e = 2 — D/2 going to 
zero, as used by Eichten and Hill in Ref.Q, then all terms such as /c^ arising from the loop 
momentum integrals vanish. In the effective theory, this implies that all contributions from 
the loop diagrams vanishes since all loop momentum integrals are proportional to it while 
in the full theory it implies that there is no any logarithmic nonlocal terms of these external 
momenta. This simplifies the matching calculations significantly. 

It is easy to see that to determine these coefficients, we need to match both the two-point 
and three-point functions. 
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A. matching two-point function 



Let's first see what we can learn from matching the two-point function. For external 
momentum p of the heavy quark near mass-shell, the general form of the QCD heavy quark 
self-energy, the inverse of the two-point function, with above infrared regulator can be writ- 
ten as: 

S(p) = A{i,-m) +2BmA, (16) 

where A is defined as 

A = (17) 

The heavy quark expansion implies that A ^ 1. Thus A, B can be expanded as power 
series of A. Up to next leading order, they can be written as: 

A = I + coifi) + C2(/i) A , 

B = cM + C4(^) A . (18) 

At one- loop level in QCD, there is only one IPI diagram contributing to the self-energy, 
as shown in Fig. 1. Carrying out a specific calculation with the above infrared regulator, 
we obtain that: 

Co(/^) = 7^as(/x)(ln-^ + 2) , ci(/i) = ^a,(/i)(ln + l) , 

C2(/i) = ^«s(/u)(ln + 2) , C4(/i) = ^as(/u)(ln-^ - 1) , (19) 

with Cf = 4/3. 




Fig. 1 Self-energy diagram on QCD side. 
The two-point Green's function reads: 

Ai4+m) - 25mA 



4m2A (^2 + AB- 52 A) 
11 



(20) 



Now we first determine the effective lagrangian up to order 1/m corrections using the 
matching condition ©. With the QCD side of © reads: 



A{2m + k-V) - 25mA 
4m2A (A2 + AB- 52 A) 



Ci + C2 + C4 



2mfc ■ 1/ 



6 



(21) 



where S = k ■ V/2m and c = 1 + Cq + Ci. 

As we argued above, on the side of the effective theory, contributions from the loop 
diagrams vanish. Thus all contributions to the right side of the matching condition (jH} arise 
from the tree diagrams. With all possible insertion of higher order terms, the right side of 
the matching condition (jHJ reads: 



1 



Zk-V 

Comparing to (ED), we see that Z = c 
the one-loop values given in (fTIJ|) . we have: 



Z, 6 



2mk ■ V 

1 + Co + Ci, and Z^ 



(22) 

(ci + C2 + Ci)/c. With 



Z 

z. 



c 



1 + -^a,(/i) 31n^ + 4 



47r 



C 



-^a.(/i)(31ni^ + l). 

2,71 ^ ^ 



(23) 



These results are in agreement with those presented in literatures. The coefficient Zm can 
only be determined by matching the 3-point function in the next subsection. 

We then use the matching condition p3|) to determine the renormalized hv^{x) field up 
to order 1/m corrections. With the two-point function given in ^B^ . the QCD side of (|TH|l 
reads: 



V G{p) F 



4m2 A (A2 + AB- 52 A) 
P_ ^P+ 



2mck ■ V 



1 + 



1 + 



C2 + C4 



• (24) 



2mA; -V \ ' c c(l + cq) 
On the effective theory side, again contributions from loop diagrams vanish while only 
tree diagrams survive. Up to next leading order correction terms five diagrams give nonzero 
contributions, as shown in Fig. 2. The Feynman rules for the operator insertions in these 
diagrams can easily be obtained from (fT^ and (fT3j) . Their contributions to the right side of 
flT^ reads: 



2m 



Zk-V 2niZ{k-Vf \2mZ 



di + d2 
2mZ 



6 



(25) 
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Fig. 2 Diagrams contributing to the matching conditions on the effective theory 
side. The circle, up and down triangles represent the operators with coefficient do, 
di, and ^2, respectively, while the solid and the blank boxes represent the insertion 
of the kinetic and {D ■ V)"^ /2m operators in the effective lagrangian. 

Comparing ()24p to ()25p . we obtain that: 

do = 1 , 

di + d2 = 1 - . (26) 

1 + Co 

We see that only the combination of the di and d2 can be determined by the matching 
through the 2-point function. To determine them separately, we need to match the 3-point 
function. 

B. matching three-point function 

In this subsection, we determine the short-distance coefficients Zm in ()14|1 and di, d2 in 
by matching the three-point function. Here the Feynamn diagrams with 3-gluon vertex 
are involved. We use the background field method^^, which can be used to simplify the 
calculations significantly. In this method, QCD Ward identity for the self-energy and the 
IPI quark-gluon vertex takes a QED-like form: 

KT>^{pi,P2) = S(p2) - (27) 

where r'^(pi,p2) is the IPI 3-vertex with external quark momenta pi, and p2, and gluon 
momentum k = p2 — pi- 

Up to next leading order correction terms, the general form of the QCD IPI vertex 
satisfying the Ward identity ()27p with quark near threshold can be written as: 

r^(pi,P2) = + ^P" + ^2P'(P-^) + S [^'^'] ' (28) 
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where 



A = 1 + Co + C2 A , 
B = ci + 2 C4 A , 

The one-loop coefficients Cq, Ci, C2, and C4 have been given in (fT^ . Thus we only need to 
evaluate the C3. For simplicity, we take the polarization vector of the gluon e satisfying: 

e ■ pi = e ■ p2 = e ■ k = 0. (30) 

Then we have: 

r(pi,P2) = T^{pi,P2)e, = A ^ + ^[fi, (31) 

In QCD, at one-loop level, there are two Feynman diagrams contributing to r'^(pi,p2) as 
shown in Fig. 3. A straightforward calculation gives that: 

cM = (2C, + C,)^{ln^^+2). (32) 




Fig. 3 Vertex diagrams on QCD side 

In the matching procedure, taking the polarization vector e can also significantly simplifies 
the calculations. In QCD, the general form of the 3-point Green's function with this vertex 
contributing to the matching conditions (jH)) and (fT^ is expressed by: 

G^(pi,P2) = G{p,)T^{puP2)G{p2). (33) 
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Now we first determine by using (jH|. The QCD side of (jHI) reads: 
P+G(pi)P(pi,P2)G(p2)P+ 



A2(^2 + 2m) - 2E2»^A2 . . 

4m2A2 (Ai + A2i?2 - 5IA2) ^ ^ 

wfiere tlie subscript 1 and 2 denote tlie momentum being pi and p2, respectively. Expanding 
it to leading order of fc, we have 

1 



I + C0 + C3) P+[^, ^]P+. (35) 



Am c^ki-V k2-V 

On the effective theory side, only the insertion of the color- magnetic dipole term gives 
non-zero contribution for the gluon polarization vector satisfying (|3(J|) . It reads: 

Zrn r . . . 

AmZki-Vk2-V +^^' ''J + ^ ) 

Comparing to dSSl), "we determine that 

= 1±^^ = 1 + . (37) 

c c 

With the coefficients given in (fT^ and (jHS), at one-loop level reads: 

47r \ m 



1 + ^ I Ca ln£^ + 2C^ + 2C^ . (38) 



We then determine di and ^2 in the renormalized hy- field fll5|l by using the matching 
condition (fT^ . The QCD side of (fT^ reads: 

p_G(pi)r(p„p2)G(p2)p+ 



4m2Ai (A? + AiBi - BfAi) V 4m 
A2(^2 + 2m) - 252mA2 



• • 4m2 A2 {Al + A2B2 - BIA2) ' ^^^^ 
The expression can be expanded as power series of /c's. Keeping only the leading corrections, 
PHl reads: 



P_ ^P ^ 



2mck2 ■ V 



1 + ^ 



2m A;2 ■ \^ 



/ 2ci - C2 - 2C3 \ Ci + C2 C4 ^ , C2 + 2C3 

- 2(1 + j ^ + 2(1^* 

1 + Co + C3 1 



4 m c2 ki-V k2-V 



F ^P+ ^]P+. (40) 
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On the effective theory side, contributions may arise the insertions of the operators both in 
hv~{x) given in (fT3j) and in the effective lagrangian given in ()14|) . With the polarization vec- 
tor e, there are 6 Feynman diagrams contributing to it as shown in Fig. 4. With appropriate 
Feynman rules, they read: 




Fig. 4 Diagrams contributing to the matching conditions on the effective theory 
side. The notations are the same as figure 2. The sohd Oval represents the insertion 
of the color-magnetic dipole operator. 



Comparing to (HH), we determine that 



di 
d2 



2Ci - C2 - 2C3 

2(1 + Co) 

C2 + 2C3 



(42) 



2(1 + Co) 

These values are consistent with ()2(i|l obtained by matching the 2-point function. 

Substituting the short-distance coefficients in (|Tn|) and (fH^ into (jl^ . we obtain the 
one- loop renormalized coefficients for di and d2: 



di{ii) 
d2{n) 



1 + 



An 



Ca - 2C^ ) In At - ICf 



"^^^^C.fln4 + 2 



47r 



(43) 
(44) 
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with Ca = 3. These are the central results of this section. These results are in disagreement 
with that obtained by Balzereit in P|. Since we use the same definition of the heavy quark 
field, both result should be equal. The results presented in this paper are derived rigorously 
using the matching conditions while Balzereit obtained indirectly from the requirement of 
the invariance of the effective lagrangian which is much more complicated. 

IV. RPI OF THE RENORMALIZED EFFECTIVE LAGRANGIAN 

In this section, I compare the renormalized transformation determined by the matching 
condition (fTSj) with those given in Q. I show that their results can easily be understood 
by constructing the effective lagrangian in an alternative way, in which an effective theory 
in four-component field is constructed first, followed by its reduction to the effective theory 
in two-component field. I then prove that the renormalized transformation determined by 
the matching conditions (|T!?)l can be written as the same form with the transformation 
(jH) and © with the covariant derivative substituted by the operator which may be called 
as the generalized covariant derivative. It means that the result presented in this paper 
is consistent with that given in j?!. Finally, I will show that the renormalized effective 
lagrangian is reparameterization invariant under the renormalized transformation. 

A. effective lagrangian in four-fermion field 

In the conventional method, a renormalized effective lagrangian is constructed by follow- 
ing steps. First a proper field to describe the low energy particles is chosen. In HQET and 
NRQCD, this effective field for describing the the heavy quark is just the two-component 
field. Then the effective lagrangian in this field is expanded as sum of local operators in 
terms of appropriate counting rules. Then the renormalized short distance coefficients of 
these local operators are determined by matching the full theory and the effective theory. 
We refer this method as matching after expansion. 

Here we introduce an alternative way to determine the renormalized effective lagrangian. 
In this method, renormalized local operators expressed in the field of the full theory are 
added to the lagrangian of the full theory by matching conditions. Then it is expanded in 
terms of the two-component field. We refer this method as matching before expansion. 
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Let's illustrate how this works in a hard cut-off regularization. 

As in the last section, we take different hard cut-off regularization energy scales Ae and 
A J in the effective theory and in the full theory, respectively. They satisfy Ag ^ m ^ Aj. 
In calculating the one-loop IPI diagrams in full QCD theory, we need to calculate the loop 
momentum integrals from zero to Af. They can be separated into integrals from to Ae 
and integrals from Ae to Af. The first part is just the same with that in the effective theory 
while the second part gives extra contributions. As argued above, the contributions from 
this region can be written as local terms of external momenta and can be expressed as 
contributions from local operators. Therefore, once those local operators are added to the 
lagrangian, the effective theory with hard cutoff Ae can produce the same result of the full 
theory with cutoff Af. This argument can easily be generalized to the case of multi-loops. 

At this stage, those local operators are written in terms of Dirac four-component field. A 
general form of the renormalized effective lagrangian density with hard cutoff Ae for heavy 
quark field can formally be expressed as: 

Ces = ^ix){ip -m)^{x) + ^{x)Oi{x)'<iJ{x) , (45) 

where = — igA^^T"- is the covariant derivative. It may be denoted as 

£efT = ^ix)Oix)^ix) , (46) 

for short by defining 0{x) = i p — m + Oi{x). 

The first term in is just the tree- level lagrangian while the second term arises from 
the renormalization with a cut-off Ae <^ m. The operators in this term are generally the 
function of the covariant derivative and the heavy quark mass. It may contain terms such 
as + m?, and QsG^^ = ^[D^, D"], which are suppressed by the off-shell momentum of the 
heavy quark or the momenta of the external gluons. They can be organized via appropriate 
power counting rules. In perturbative calculations, the loop momentum integral is from zero 
to Ae. In this region, both the external and the loop momenta are smaller than m, hence 
the 1/m expansion is allowed and the quark mass dependence is extracted explicitly. Thus 
the energy scale m is no longer involved in the effective theory. I emphasis here that the 
effective lagrangian density in this form is independent of the velocity parameter V. Thus 
it automatically satisfies the RPI. 

Higher dimensional operators appear in Oi{x). It implies that power divergences arise in 
the loop momentum integrals. In the full theory the power divergences cancel when both the 
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contributions from quark and antiquark are included. However, in the effective tlieory when 
we impose a hard cutoff Ag <^ m on the loop momentum integrals, the contributions from 
antiquark are excluded so that the power divergences do not cancel. Nevertheless, those 
power divergences are artificial since the power divergences from the diagram calculations 
just cancel those from the short distance coefficients. 

At one loop and leading order of 1/m, the most general form of the four- component 
effective lagrangian is 



eflf 



'1 + Co)^(x) (i/) -m)^(x) - —^(x){D^^ir?)^{x) 

2m 



«C2 ^ , . 



C3 



Am 



i p - m){D^ + m^) + {D^ + m"^ ) {i p - m 



5m'^ 



(47) 



Calculating the the IPI diagrams shown in Fig. 1 and Fig 3 using this effective lagrangian 
and full QCD, we see that these coefficients are just those Ci — C4 given in (|T^ . and ((221) • 



B. effective lagrangian in two-component field 

Now let's reduce (I47|l to the effective lagrangian in two-component field. The equation 
of motion now reads: 

P-0{x) {hv+{x) + hv-{x)) = 0, (48) 

where 0{x) is the 0{x) whose covariant derivative iD is replaced by iD + mV due to the 
phase factor in the field redefinition. It can be regarded as the renormalized equation of 
motion. 

From Eq. ()48p . we can express hv-{ clS db function of hv+{x) formally as: 

hv-ix) = ^ -nr ^ T/ V n( P- /P + Oi{x)) hy^{x) , (49) 
2m + iD[x)-V — r_C'i(x)r_ ^ ^ 

where Oi{x) is the Oi{x) whose covariant derivative iD is substituted by iD + mV . This 

modifies the tree level expression Once the form of Oi{x) is given, the right hand side 

of Eq.(jini) can be expanded as power series of 1/m. With Oi(x) given in (jUj), up to order 

as and 1/m^, hv-.{x) reads: 



hv-{x) 



2m ^ 4m2 I 2(1 + Cq) / 



C2 + 2C3 



8(1 + co)m2 
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pD ■ V hv+{x) . (50) 



Comparing this with ()15p. ()42|) . we see that they are in agreement. 

Finally, with equation of motion (jIH|) . the effective lagrangian (jlUj) is reduced to: 

£eff = hv+{x)0{x) (^hv+{x) + hv-{x)'j 

= [hv+{x) + hv-{x))0{x)hv+{x) . (51) 

This is just the two-component effective lagrangian. It can be expanded as power series of 
1/m. In this way, the four component effective lagrangian is reduced to the two component 
effective lagrangian. Up to order and 1/m correction, it reduced to the effective lagrangian 
fll4|) . Therefore, the effective lagrangian obtained by these two different approaches are just 
the same. 

C. comparison with the previous work 

In Sec. II, the matching condition for determining the renormalized hy_{x) field is given 
by Eq. (fT!^ . In the last subsection, the hv-{x) field has been obtained by equation of motion. 
Its expression is given by Eq. (j49j) . In this subsection, I will show that the small component 
fields obtained by these two different methods are identical. They uniquely determine the 
renormalized transformation of the heavy quark field against the infinitesimal variation of 
the velocity parameter V. Adding both sides of (fT^ and (fTSj) together, we have 

{0\T[hv+{x) + h'y_{x))hvAy)\0)'' = (0|TM/(a;)^(y)|0)«P+. (52) 

where (0|T\I'(x) 'i'(?/)|0)^ is a full propagator under arbitrary external field B^{x). It is 
satisfied order by order in as- Suppose we calculate the left hand side at tree level with the 
renormalized effective lagrangian. To validate this equation, the right hand side then should 
also be calculated to the tree level with containing contributions of the loop momentum 
integrals from Ag to Aj. This can be calculated by the renormalized four- component effective 
lagrangian to tree level. Thus it satisfies the following equation 

0{x)G{x,y-B) = i6\x~y). (53) 

Acting an operator P_0{x) on the left hand side and P^O{x) on the right hand side of 
the right hand side vanishes immediately since P_ ■ P+ = 0. Since we only calculate them 
at tree level, the operator 0{x) can be moved within the bracket: 

P.{0\T[O{x){hv+{x) + h'y_{x))hv+{y)]\0)'' = 0. (54) 
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Since the argument y in hv+{y) is arbitrary and this correlation function contains interaction 
with arbitrary background gluon field, the unique solution of this equation is 

V.O{x)(^hv+{x) + h'y_{x)) = 0. (55) 

This is just identical to ()48|) if hy_{x) is the same with hv-{x). This implies that the 
renormalized hy__{x) determined by the matching condition ()13|) is identical to that from 
the equation of motion (j48|l . 



D. RPI of the renormalized effective lagrangian 

In this subsection, let's prove that the renormalized effective lagrangian ()5H1 is invariant 
under the transformation (j3]) or (jllll with renormalized small component field. 

It follows that from an infinitesimal transformation of the effective lagrangian 1)511) 

A£eff = AhY^{x) 0{x) hv{x) + hv+{x) 0{x) Ahv{x) 

.A^— - — 

= hvix)^-0{x)hvix) + hv+ix)0{x)Ahv{x) . (56) 

We have use a shorthand notation hy{x) = hv+{x) + hy_{x). It is emphasized here that the 
operator 0{x) which is from the four component effective field theory is invariant against 
the variation of the velocity V. Any change arising from the phase factor in the definition 
of the effective field has been omitted simply because it is trivial under the transformation. 
Imposing an infinitesimal transformation on the equation of motion (^Hj), we obtain that 



- ^0{x)hvix) + F-0{x)Ahvix) = 0. (57) 
With it, ()5(jj) can be rewritten as 

A^cff = hvix)0{x)Ahvix). (58) 
Notice that P+ hv~{x) = 0. Imposing an infinitesimal on it, we immediately have 

P+A/iy_(x) = -^hv-{x). (59) 
Adding it together with P^Ahv+{x) = AV'/2 hv-{x), we have 

F+Ahv{x) = 0. , (60) 
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With it, is reduced to 

A£eff = hv{x)0{x)F^Ahv{x) . (61) 

It follows that ACcs = from the equation of motion hv{x) 0{x) P_ = 0. Thus we have 
shown that the renormalized effective lagrangian (j3T|) is invariant against the variation of 
the velocity parameter V under the infinitesimal transformation (ji} with the renormalized 
small component field. 



V. CONCLUSION 



The RPI is an important theoretical issue in the heavy quark effective theory and the 
NRQCD effective theory. It is required by the consistence of the effective theory. It also 
leads to interesting applications 10| 11| 3|- The transformation of heavy quark field under 
the variation of the velocity parameter V proposed by Chenj^ with tree level expression of 
the small-component field keeps the tree level effective theory invariant. However, at loop 
level, the transformation needs to be renormalized in the renormalized effective theory. In 
this paper, I show that the renormalized transformation of the heavy quark keeps the same 
form as Chen's transform while the small component field needs to be renormalized. I derive 
the matching conditions for determining the renormalized transformation by imposing an 
infinitesimal transformation on the relations between the Green's functions in the full QCD 
and those in the effective theory. These matching conditions are essential for studying the 
renormalization issue in RPI. As an application of these matching conditions, I determine 
the renormalized transformation up to order I also show that the previous result 



in |3| can be understood clearly by building the effective theory in an alternative way, in 
which the a renormalized effective lagrangian in Dirac four-component field is constructed 
first, followed by its reduction to the two-component effective lagrangian. The renormalized 
small component field is then obtained by the equation of motion. The four-component 
effective lagrangian automatically satisfies RPI. Thus RPI cannot give any constraints on 
any operators in it. When it is reduced to the two-component effective theory, the same 
operator with certain coefficient may appear in different terms. The RPI can be used to 
connect those terms. I also show that the obtained renormalized small component fields 
by these two methods turn out to be equivalent while the matching conditions provide a 
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systematic way to determine the renormalized transformation to any desired order in 1/m 
and as expansions.. 

APPENDIX A: GENERATING FUNCTIONAL OF GREEN'S FUNCTIONS 

In this appendix, we derive the relations between the Green's functions in QCD full 
theory and that in the effective theory using generating functional method. It is similar 



with that given in Q and 0. We use the background field method 
interactions, which preserves explicitly the gauge covariant. 
In QCD full theory, the generating functional reads 



14 



3 



for gluon field 



Z[ri,ri,J,B] = J 5[ri,r],A]expi J d^x{lQ{x) + Ig{x)) 



(Al) 



where rj, ri,J are the external sources for heavy quark, antiquark and gluon field, B is the 
background gluon field, Ig is given by 



Ig(x) = --F^^^F"'"' - —{Cy + Indet 



with 



(A2) 



F;, = d,{A + B)l -dM + B)l + gr'^A + B)%A + B) 



(A3) 



= d,Al + gr^'BlA'^^ , 



(A4) 



being the gauge-fixing term. If satisfies the following relation 



5W 

Jb^. 



+ d'y 



SWSJtiy) 



6Jt 6B- 



- 7" 



(A5) 



with W[ri, T], J, B] = —i In Z[ri, r], J, B], W[ri, t], J, B] is just the effective action regarding to 
the gluon field B with gauge-fixing term 



= dM-B)l + gr''BlAl, 



(A6) 



and In reads 



Iq{x) = '^{x){i p - m)'^{x) +r]{x)'i!{x) + '^{x)'r]{x) 



(A7) 
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The quark field can be integrated out formally and then we have 

Z[ri, f], J,B] = J 5[A] det[i p~m\expi J d^x{lQ + Ig) , (A8) 

where I'q remains the same form as Iq. But the quark field now is related to the external 
source r]{x) by the following equation of motions: 

{i P - m)\l>(x) = -ri{x) . (A9) 

The generating functional of the effective theory is similar with that of the full theory 
except the heavy quark action. The effective lagrangian is substituted by ()46|1 . In the 
external source term of the heavy quark only the large component effective field defined in 
(HI) couples to the external source. The action of the heavy quark is given by: 

= ■^(x)0(x)^(x) + f]{x)V+hv+{x) + /iy+(x)P+r/(x) . (AlO) 

Similarly, integrating out the heavy quark field, the generating functional takes the same 
form as ()A8|) with the effective action of the heavy section is substituted by 

I'^^{x) = hv{x)0{x)hv{x) + f]{x)F+hv+{x) + hv+{x)F+v{x) , (All) 

with hv{x) = hv+{x) + hv-{x). 

The quark field now is related to the external source ri{x) by the following equation of 
motions: 

0{x)hv+{x) = -F+vix) . (A12) 

Multiplying P_ on both sides, the right hand vanishes and we obtain the renormalized 
equation of motion: 

P_ 0{x) [hv+ix) + hv-ix)) = . (A13) 

This is just the equation (jlHj) . The renormalized hv-{x) can be related to hv+{x) by 
With the equation of motion ()A13|1 , (jAlljl can be simplified as 

Iq = hv+ix)0{x)hv{x) + f]{x)P+hv+{x) + hv+{x)P+r]{x) . (A14) 

This gives the effective lagrangian density (|5ip. 
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The quark determinant in ()A8|) is responsible for the contributions of the heavy quark 
loop. It is the same in the full theory and in the effective theory and is suppressed by 1/m? 
at least. Thus we may ignore it. 

The full quark propagator with background field B^{x) is gained by differentiating over 
the external sources. 

Gix, y; B) = PFfa 7^, J, B) . (A15) 

dr]{x)dri{y) 

If the hard cutoff energy scale is set to Af, the same with that in the QCD full theory, 
the 0{x) is then to be z p — just the same with that in the full theory. The effective 
lagrangian is just the nonlocal form 0. In this case the only difference of the effective 
theory and the full theory is the external source term. One immediately gains the relations 
between the Green's functions of the full theory and those in the effective theory (jH)). This 
relation ensures that the nonlocal effective theory is equivalent to the QCD full theory. The 
local effective theory with the hard cutoff regularization scale Ag is equivalent to that the 
nonlocal effective theory with hard cutoff Aj. This ensures the validation of the relation (jHJ. 
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